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ANALYTICAL PROCEDURE FOR EVALUATING 
SPECKLE -EFFECT INSTRUMENTATION 
by Arthur J. Decker 
Lewis Research Center 

SUMMARY 

A general analysis suitable for developing speckle-effect instrvunents and a simpli- 
fied analysis suitable for evaluating laser -speckle instrumentation are presented. The 
simplified analysis is svunmarized as a list of equations. Existing speckle instrumenta- 
tion and techniques can be understood by an application of the principles of geometrical 
optics and the orderly use of these equations. To demonstrate the use of the simplified 
analysis, several sample applications are described. Two forms of speckle interfer- 
ometry are analyzed: one form is used to measure strain, and the other form to display 
the nodes of vibration of an object vibrating in one of its eigenmodes. A simple speckle 
viewer and the propagation of spatially coherent light from a rough surface are also eval- 
uated, The properties often stated for speckle and the assumptions that warrant these 
statements are reviewed. 


INTRODUCTION 

In the future, the speckle effect probably will be applied in aeronautics research and 
development as a valuable supplement of holography. "Speckle effect" refers to the 
granular appearance of a rough surface illuminated with spatially coherent light. The 
large variation in intensity observed from speckle to speckle is caused by the random in- 
terference of the coherent light reflected from the rough surface. The apparent typical 
size of the speckles is approximately equal to that of a resolution element of the optical 
system used to review the rough surface. Speckle can be used in the measurement of 
surface strain (ref, 1), in the measurement of surface roughness (ref. 2), and in the 
measurement of the amplitude distribution and node locations on a vibrating object 
(ref. 3). In the measurement of surface strain, the speckle effect can be used to deter- 
mine strain in the plane normal to the viewing direction. Holographic techniques can 



yield strain approximately along the viewing direction. Hence, the speckle techniques 
and the holographic techniques are definitely complementary for the measurement of 
strain. In determining the node locations on a vibrating object, vising a speckle interfer- 
ometer eliminates the need for the intermediate photographic step required in real-time 
holographic interferometry. Furthermore, as reviewed in this report, speckle interfer- 
ometry can be used to yield time-average fringes for measuring the amplitude distribu- 
tion of a vibrating object. These fringes appear within the integration time of an eye or 
a television vidicon and do not require the photographic step required for time-average 
holography. However, the visibility of the fringes formed in speckle interferometry is 
not nearly as good as the visibility of the fringes formed in time -average holography. 

A guide is needed for analyzing the speckle phenomenon that can be used by those 
who wish to design, apply, or simply purchase speckle instrumentation for use with a 
variety of surfaces. Unfortvmately, simple mathematical descriptions of the speckle phe- 
nomenon and its applications are not as easily generated as are descriptions of the holo- 
graphic process. Indeed, a substantial part of speckle -effect research is devoted to de- 
termining the field and intensity statistics of spatially coherent light reflected from 
rough surfaces described with different roughness models (refs. 4 to 6). 

This report evolves a guide for a general analysis of speckle -effect instrumentation 
based on standard coherent optics techniques (ref. 7). Then, the analysis is cast into a 
form that is not sensitive to surface properties and that yields some of the conclusions 
frequently quoted for speckle. These conclusions include a speckle size that is deter- 
mined by the resolution of the optical system and a distribution of intensity that is a neg- 
ative exponential. Where simplifyii^ assumptions are required, they are carefully 
stated. Thus, the more general analysis can be referenced when the assumptions are 
not valid. Three examples of the application of the analysis are presented. First, the 
results for simple propagation of radiation between a rough surface and a detecting plane 
originally derived by Goodman (ref. 8) are derived. Next, the analysis is applied to a 
simple speckle viewer. Finally, a speckle interferometer of the type described by 
Stetson is analyzed (ref. 3). The intention is that this report can be used by those who 
wish merely to evaluate speckle instruments and techniques. However, more complex 
effects of instrvunent components or deviations of the speckle effect from the norm can 
also be analyzed. 


GENERAL ANALYSIS OF A SPECKLE INSTRUMENT 

The speckle instrument is treated as a linear system that transforms the light orig- 
inating at a rough reflecting svmface to a distribution of electric field in a detectii^ plane. 
The light is assumed to be spatially and temporally coherent. The rectangular compo- 
nents of the electric field are assumed to be uncoupled: They satisfy separate wave 
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equations and separate transformation equations. If E^(x,y) and Ey(x, y) are the 
phasor amplitudes of the x- and y -components immediately after reflection from the 
rough surface and E^{x,^ and E^(x^y) are the phasor amplitudes in the detecting plane, 
the linear transformation relating the fields is (ref. 7) 


E^(x,y) 



x,y)Ej^(x,y)dx dy 


-OO 



OO 

■// 


hy^,y;x,y)Ey(x,y)dx dy 






( 1 ) 


where h^ and hy are the impulse response functions for the x- and y-polarizations of 
the field. For complex optical systems, h^ and hy can be constructed from a sequence 
of linear transformations. These transformations are easily written for thin lenses with 
aberrations and finite dimensions and for simple propagation of light. Polarizations are 
specifically noted since the impulse response function may depend upon polarization. 
Frequently, only a single polarization is admitted (ref. 3). Then, either h = 0 or 

The field reflected from the rough sxirface is treated as a stochastic variable. Po- 
larization, magnitude, and phase at a point x, y are random quantities characterized by 

a probability density function f(E E ). Only the statistics of the field are predictable. 

X y 

As stated in the INTRODUCTION, the specification of the probability density function f 
is the subject of continuing investigation. 

The goal of this report is to develop an analysis that will apply to a variety of den- 
sity functions. Typically, (E^) = 0 and (E~> = 0, where () denotes an ensemble aver- 
age. In defining () , one imagines a set of an infinite number of macroscopically iden- 
tical surfaces where the local microstructure varies from one surface to another. Then, 
the symbol () refers to an arithmetic mean over all members of the set. The first- 
order statistics of the field are of no value in understanding speckle effect. However, 
the second-order statistics in the field and in the intensity are useful (ref. 8). In par- 
ticular, the intensity correlation or the magnitude of the field correlation tends to be high 
within a speckle. Using the notation of reference 8, the correlations in the detecting 
plane are 
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OO 






xh^(xj,yi;x2,y2)<E*(Xj,y^)E^(x2,y2))dXj dy^ dj^ dy^ 


> ( 2 ) 


OO 

r^(^2’y2’^i'yp = ^ jfj f 


><hy^l»yi;x2,y2)<Ey(Xj,yj)Ey(x2,y2)>dx^ dy^ dx2 dy2 




Goodman in reference 8 calls these correlations the coherence functions. When nearly 
monochromatic radiation is used, these functions are a useful extension of the concept of 
coherence from the temporal domain to the spatial domain (ref. 9). Coherence theory 
has traditionally considered time averages or time as a parameter in evaluating an en- 
semble average (refs. 9 and 10). However, the nature of speckle instrumentation and 
the availability of highly monochromatic laser radiation make these spatial coherence 
functions more useful. Ensemble averages are evaluated with position as a parameter. 

The actual relation between and Pj;; depends on the details of polarization. The 
ensemble average <Ej^(Xj,yj)E^(x2,y2)> depends on a density function fJE^^ ’^x« 

3k — ' 1 "2/ 

similarly, the ensemble average (Ey(xj^,yj^)Ey(x2,y2)) depends on a density function 
fy^y where 1 and 2 in the svibscripts refer to the positions (xj^,y^) and (x2,y2). 

For completely random polarization, does not depend on the y -polarization and 
does not depend on the x-polarization. Two completely independent speckle patterns 
ensue and are superimposed. 

For a completely deterministic polarization. 




The function f 


f (E ,E ) is used to express the probability density function for an x- 
^1 ^2/ 

polarized field E^ at position (xj,y^) and a x-polarized field E^ at position (x2,y2) 
1 2 
in terms of a function of the corresponding y-polarizations. Similarly, the function 


4 



f,iE E ) is used to express the probability density fmction for a y -polarized field 

Xi Xg/ , _ 

E at position (xj,yj) and for a y -polarized field E at position (Xg,y 2 ) in terms of 

i a 

the corresponding x-polarizations. 

The quantity <E*(xj,y^)E^(x 2 ,y 2 )> is proportional to <E*(xj,yj^)Ey(x 2 ,y 2 )> • If h^ 

is proportional to hy (simple attenuation and phase shift of one polarization relative to 

the other), r~ is proportional to 1^. For circular polarization and li = h = r~. 

In any case, for a completely deterministic polarization, a single speckle pattern is 

created and either or r~ svtffices for evaluating speckle size. 

X y 

The analysis that immediately follows assumes that the polarization is deterministic 
or that the instrument has a polarizer that sets hy = 0. A polarizer is easily added to a 
speckle viewer and is required for the speckle interferometer discussed imder EXAM- 
PLES (ref. 3). The case of random polarization is discussed in the section Intensity 
Probability Density and Overlapping Speckle Patterns. 


Speckle Instrument: Single Polarization or Light Deterministically Polarized 
The subscript x will be dropped in referring to E and h The subscript x will 

X X 

be dropped in referring to r^. The field correlation or spatial coherence function in the 
detecting plane is then 


CO . 

r(x2,y2’^i»yi^ = //// h*(x2,y2’*l’yi^^^pyi’*2’^2^ 


x(E*(xj,yj)E(x2,y2))dXj dyj dx2 dyg (3) 

The entire space between the rough surface and the detecting plane is called the optical 
system. This system is divided into three parts. The key part is a physical aperture 
stop that defines both the exit and entrance pupils. This stop may be, for example, an 
iris of adjustable size, a lens, a roxigh surface of finite extent, or a detector of finite 
extent at the detecting plane. The stop is defined by the equations 


P(x„,y„) = 1: when (x„,y„) is a point in the aperture 

P(Xa, y^) = 0 otherwise 

The optical system between the physical stop and the rough sxmface is described by the 
impulse response function 


(4) 
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y) 


where f refers to the front section. The optical system between the physical stop and 
the detecting surface is described by the impulse response function 



where b refers to the back section. The impulse response function for the entire sys- 
tem can be shown to be 


h(x,y;x,y) 




X, y; x^, ya)P(x^, Y 


(5) 


The quantity hj^ • P(x^)Ya) primarily responsible for typical speckle size. Be- 

fore writing r again, make the following changes of variable. 


II 

Yj = y 

X2 = X + Ax 

Y2 = y + ^y 


ya = ^a 


Ya 2 =ya + ^ 

Xj = X 

Yi = y 

r>j r«w> 

Xg = X + Ax 

yg = y + Ay 


a 


These variable changes and equation (5) are substituted in equation (3) to give 


+ Ax,y + Ay;x,y) = <E~(x + Ax,y + Ay)E~(x,y)> 

'//////// 




X hj,(x, y;x^ + Ax^, y^ + Aya)P(Xa, ya>P(x^ + Ax^, ya + ^ya> 

X hf (x^, y^; x, y)hj(x^ + Ax^, y^ + Ay^; x + Ax, y + Ay) 

X <E*(x,y)E(x + Ax,y + Ay))|dx dx^ d Ax d Ax^ dy dy^ d Ay d Ay^ (5) 
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The pairs of variables x and y, x and y, and x and y and their subscripted and in- 
cremented forms appear symmetrically in many of the equations that follow. To shorten 
the equations, only a single variable is written to stand for both variables in the pair. 
Define the correlation in the aperture stop as 


//(-•' 




AXg^; X + Ax)( E 


^ 


(x)E(x + Ax)) I dx d Ax (7) 


where x only is written to stand for x,y. When equation (7) is substituted into equa- 
tion (6), 


OO 

+ Ax;x) = AS; + Ax^) <E^(x^ + Ax^)E*(x^)> dx^ d Ax. 


( 8 ) 


Again dx„ d Ax„ is taken to mean dx„ dy„ d Ax„ d Ay„. With a finite aperture, the 
limits of Ax^ depend on x^. 

The expression for the spatial coherence function is seen to be quite general. The 
coherence function in the detecting plane is a transformation of the coherence function in 
the apertrare. The transformation is four -dimensional depending on x„, y_. Ax , and 
Ay^. Two assumptions are required if a simple interpretation of the speckle character- 
istic size is to be made. The first assumption is that 

is at least locally wide -sense stationary. That is, the normalized aperture correlation 
depends only on Ax^ but not on x^ at least over the dimensions of the aperture. The 
second assvunption is that h^^ + Ax;x^) x h^{x;x^ + Ax^^) can be broken approximately 
into three factors: 

h^(x + Ax; x^)h^^; x^ + Ax J « D(x; Ax)Tp(AXg^; x)T Jx^; Ax) (9a) 


The reason for this factorization is made apparent later. This assumption may occa- 
sionally be a poor assumption; however, it simplifies the resiilts and facilitates the un- 
derstanding of speckle phenomena. The coherence function becomes 
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(9b) 


r(x + Ax;x) = / — _ 

J_ 


Tj.(Ax^;x)d AXjj 


oo 

^_nr» 


In equation (9b), r is expressed as a product of a two-dimensional transformation of the 
normalized correlation or the spatial coherence fimction of the aperture field, a two- 
dimensional transformation of the aperture intensity, and a factor that depends on x and 
Ax. This factor D^, Ax) usually consists of phase factors if the Fresnel diffraction ap- 
proximation is made. 

If the back optical system consists only of Fraunhofer diffraction from the aperture, 
the decomposition into three factors is exact. The coherence function F is then the 
product of the spatial power spectral density of the normalized correlation of the aper- 
ture field and the Fourier transform of the intensity in the aperture. 

A normalized coherence function is more useful for evaluating speckle properties. 

It is assumed that 


/ 


<Ea(Xa + ^a^^a^^a^) 

<E:(^a)Ea(^a» 


Tp(Ax^;x)d Ax^ 


is approximately constant over the detecting surface. This assumption is generally 
equivalent to assuming that 




is a narrow function of Ax . Normalizing with respect to ' Ax = 0, 
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y(x + Ax;x) = 


r(x + Ax;x) 
r(S;x) 


_ D^; Ax) 
D(x; 0) 








a 


(10) 


Equation (10) is a primary result of this analysis. The equation expresses the degree of 
spatial coherence of the field in the detecting plane. In the case of constant intensity in 
the aperture, the degree of spatial coherence depends only on a transformation 
Tg^(Xg^; Ax) of the pupil fimction P(Xg^, y^^) defining the physical aperture stop. A number 
of assumptions have been required to reach this point in the report. The correctness of 
the assumptions must be determined on a case -by -case basis. 

A number of properties are to be noted, 

y(x + Ax;x)D(x;0) 

, D(x; Ax) 


or |y| depends only on Ax and is, therefore, stationary in the wide sense (ref. 11). If 
D(x; Ax) is a phase factor involving x and Ax, the following properties will characterize 

|rh 


y < 1 


|y(x + Ax;x)| = |y(x-Ax;x)| = |y(Ax)| 
|y(0)| > [y(Ax)| 

|y(0)| =1 






( 11 ) 


Intensity Correlations and Speckle Area 

The quantity actually detected in the detecting plane is the intensity rather than the 
field. Speckles actually detected are regions of high correlation between intensity fluc- 
tuations rather than regions of high correlation between fields. That is, if I^) - (I(x)) 
is an intensity fluctuation, that quantity shovild be consistently positive (bright speckle) or 
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consistently negative (dark speckle) within a speckle. The field has been used since the 
optical system is represented as a linear transformation on the field rather than on the 
intensity. Fortunately, the typical extent of the region in which fields are well corre- 
lated is identical with the typical extent of the region in which intensity fluctuations are 
well correlated (refs. 8 and 10). The correlation between intensity fluctviations is 

< [l^ + Ax) - + Ax))] [l^) - <I^)>]) = (I^ + Ax)I^)> - <I^)> (I^ + Ax)> 

A normalized intensity fluctuation correlation is 



+ Ax)I^)) - + Ax)) (I(x)> 

<l2(S)) -(I(x))2 


(12a) 


provided that the statistics of the field are Gaussian (refs. 8 to 10). The statistics are 
discussed later. 

A typical speckle dimension at ^,y) should be approximately the interval Ax or A^ 
over which |y| is significant. However, jyj generally will be nonzero over the entire 
detecting surface. For example, a square rough surface of side I viewed at a distance , 
d in the Fraunhofer diffraction region results in 

r*-j rsj 

I Ax . I Av 
sine — sine 

Xd Xd 



where sine x = sin ttx/ttx. This result is derived later. 

One estimate of speckle size is that it is approximately given by the first zero of the 

sine function: Ax « Xd/Z . This estimate is confirmed if a definition of speckle area is 
s 

used that is derived from a physical argument given in reference 8. The physical argu- 
ment assumes Gaussian statistics. Speckle area is defined as 


Ag=A yyd Ax d Ay|y(x + Ax,y + Ay;x,y)|‘ 


(12b) 


where A^ refers to the area of a typical speckle and A^ is the area of the detector at 
the detecting surface. For purposes of instrument evaluation and this report, A^ is as- 
sumed to be large when compared with A . Equation (10) is used for the evaluation of 

O 

r- 
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Then, 


A 


s 


= 4 


oo 



(13) 


where |y(Ax, Ay)|^ is simply used as a weighting function. If |r| =1, the whole detec- 
tor area is covered with a single speckle (no granular appearance). If |y| « 0 except 

very near Ax = 0 and Ay = 0 (low spatial coherence), A will be very small and gran- 

s 

ulation will not be resolved. Intermediate cases will give an A_ smaller than the de- 
tector area but large enough to be resolved. 

Equations (10) and (13) are the most important for instrument evaluation. 


Intensity Probability Density and Overlapping Speckle Patterns 

The intensity at the detecting plane is characterized by a probability density function. 
Although there is a high -intensity correlation within a speckle, the intensity from speckle 
to speckle will vary at random. Only a large variance in intensity yields the visual high- 
contrast speckle pattern. The exact probability density function will depend on detailed 
surface statistics. Application of the Central Limit Theorem is assumed in order to 
avoid having to know the detailed surface statistics. 

Referring to equation (1), an impulse response function of a typical optical system 
superimposes the fields from many correlation elements on the rough siirface. Again, a 
correlation element is a region with respect to Ax, Ay over which ( E (x + Ax, y + Ay) 
E^(x,y)) has a significant value. It can be shown that a minimum correlation interval is 
1 wavelength of the light used. Outside a correlation region, (Ej^(x + Ax,y + Ay) 

E*(x,y)> *=« 0. This report goes one step further and assumes that the fields are effec- 
tively statistically independent outside a correlation interval. Also, <E^(x,y)) =0. In 
general, the integral in equation (1) is divided into a sum of integrals, each integral eval- 
xiated over a correlation region of the rough surface. These integrals are taken to be 
zero mean and statistically independent. The Central Limit Theorem (ref. 11) can then 
be invoked to show that in the limit of an increasing number of correlation elements, the 
integral in equation (1) will have a normal distribution. An exact proof of this conclusion 
for a random phase diffuser is given in reference 12. If (I) is the mean intensity in. the 
detecting plane, the intensity will have the probability density function (ref. 8) 


(14) 


P(I) = Jl 

(I> 
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The variance of this distribution is (I> . Hence, the standard deviation equals the mean 
intensity and the high speckle contrast is easily understood. Equation (14) is included 
with equations (10) and (13) for the evaluation of an instrument. 

Various applications and manifestations of the speckle effect involve a superposition 
of two speckle patterns. As stated previously, if a rough surface that reflects randomly 
polarized light is viewed with a speckle instrument that admits both polarizations, two 
independent speckle patterns are superimposed. If two identical rough-surface images 
are viewed superimposed and with slight relative displacement, this displacement or 
strain can be measured by using the right instrument (ref. 1). 

First, consider the case of randomly polarized light. Write the intensity fluctua- 
tions as 


Aly=ly(j)-(ly®> 

where refers to the intensity in the x-polarized pattern and ly refers to the inten- 
sity in the y-polarized pattern. From equation (12a) and equation (14), 

+ Ax)Al^(x)> = <Ix(x)>^|yxl^ 

<Aly^ + AS)Aly(?)> =<ly(x))^|yy|2 

The fluctuation of the total intensity is AI = AI^^ + Aly = I^(x) + Iy(y) - (^(x) + ly^)) . 
Since AI^^ and Aly are independent, zero-mean, random variables, the correlation of 
the fluctuation of the total intensity is the sum of the correlations of Al^ and Alyi 

<AI^ + Ax)Al^)> =<lx^))^|rxl^ + <ly(x)>^|yyP 

Similarly, (I^^)> - <I^))^ = (Ijj(x^)^ + <Iy^)>^. Hence, the normalized correlation of 
the total fluctuation is 

Equation (15) is used in equation (13) to determine typical speckle size when the 
polarization is random. K the average intensity in the two patterns is the same and if 
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the optical system affects both polarizations equally, [y^|^ ~ ~ 

speckle size is the same as for deterministically polarized light. 

Next, consider the two identical but slightly displaced speckle pattern images. This 
case occurs with one of the more popular forms of speckle pattern interferometry 
(ref. 1) and is included for completeness. The intensity in one of the patterns is desig- 
nated I^) and the intensity in the other pattern is designated I(x + 5), where 6 is the 
displacement. The optical process, which may include photographic steps, is chosen to 
superimpose the two speckle patterns so as to produce an output field proportional to 
l(x + 6) + I^). The Fraunhofer diffraction pattern of this field is then viewed at some , 
distance d. The field viewed is 


Ep,^) cc 


/ oo ^ 

[i® +i6;+ a]e-i(2™A<»dx 


. oc 


f ^ ^i(2ra/Xd)6] 


(16) 


The intensity detected in the Fraunhofer diffraction pattern is 


1 -00 ^ 
y I(x)e-j(2^^x/Xd)^ 


[2 + 2 cos 
I Xd ; 


(17) 


The viewed intensity consists of a factor containing the Fourier transform of the speckle 
pattern intensity multiplied by a sinusoidal factor whose frequency is proportional to the 
displacement 5. The only effect of the stochastic variable I^) is to create a field with 
a broadband Fourier transform. With a broadband transform, there will be illumination 
over many cycles of [^2 + 2 cos (2ff6/Xd)xJ in the Fraunhofer diffraction pattern. Thus, 
an accurate measurement of displacement 6 is possible. The transform 




will change from one sample function I^) to another. Thus, the Fraimhofer diffraction 
pattern Ip(x) will change from one sample of the rough surface to another. In order to 
establish a unique transform, (Ip(x)> is usually dealt with. 
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If the sensitive area of the detector at the detecting plane is large, 


. a; 

<1^ + Ax)I^))e"^^^^ Ax/xd)^ 

nr» 


Ax (2 + 2 cos ^IL£x 
X d 


( 18 ) 


From this resxilt comes the confusing statement that this type of speckle interferometry 
displays the Fourier transform of the autocorrelation function of a speckle pattern. The 
physically useful fact is that the rapidly varying intensity in a speckle pattern assures , 
uniform illiimination in the Fraxmhofer diffraction pattern for ease of measuring 
displacement. 

The general analysis is concluded at this point. 


SUMMARY OF EQUATIONS 


Before the analysis is applied to examples, the essential equations are summarized; 

(1) Decomposition of the optical system into a front part designated f, a physical 
aperture stop ^nd a back part designated b; 


00 

h(S,y;x,y) = j J h^,(x,y;x^,y^)P(x^,yJhj(x^,y^;x,y)dx^ dy^ 


(5) 


(2) Approximate factorization of h^^ + Ax;x^)h^^;x^ + Ax^) into three factors: 

hjl^ + Ax;x^)h^^;x^ + Ax^) « D(x, Ax)T^(Ax^;x)T^(x^; (9a) 

where x is a shorthand for ^,y) and x_ is a shorthand for (x„,y„) 

' H- cL' ci 

(3) Field correlation or coherence function in the detecting plane expressed in terms 
of the field correlation in the physical aperture: 


/ OO 

<E?^a)^a(^a» 


Tj,(AXg^;x)d Ax^ 
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(4) Spatial degree of coherence in the detecting plane: 


y(S + Ax;x) = 

D(x; 0) 




L 


<E:(Xa®a<='a»Ta(Xa;0)dX3 


( 10 ) 


(5) Intensity fluctuation correlation for Gaussian statistics: 

|y|2 _ <I^ + Ax)i^)> - (1^ + Ax))(I^)> 

which is equation (12a) rewritten with a slight change since (I^(x)) = 2(I(x))^. 

(6) Typical area of a speckle on a large detector; 


00 

■•// 


d Ax d Ay j y(Ax, Ay) | ^ 


(13) 


(7) Probability density of intensity when central limit theorem applies: 


P(I) = J_ e"^Al> 

<I> 


(14) 


(8) Normalized correlation of intensity fluctuations for randomly polarized light: 




2 ^"x 




(15) 


(9) Ensemble average of intensity in the Fraunhofer diffraction pattern of overlapping, 
identical, but displaced speckle patterns: 


/.CO ^ / \ 

(Ip(S> oc y (I^ + Ax)I( 5^)) AxAd)^j ^3 + 2 cos ^ xj (18) 
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EXAMPLES 


Propagation of Light between a Rectangular Rough Surface 
and an Infinite Detecting Plane 

Figure 1 shows a rectangular surface of size L by L separated by distance D 

X y 

from a detecting svtrface of infinite extent. In this case, the aperture stop is the object 
itself and is specified by the function 


— 

^ ^ L L 

X y 

Rect X = 1 when |x[ <1/2 and rect x = 0 when |x| > 1/2. The factors in equations (5) 
are easy to write down. 

If Fresnel diffraction is used to describe propagation (ref. 7) 

hb<s.y;='a-ya> H^a -^^1} 

A phase factor has been dropped. In preparation for the factorization of eqxiation (9a), 


h^^ + Ax;x^)h|3^;x^ + Ax^) = 



X exp 1^- ^ (x Ax^ + y Ay^)j exp (x^ Ax + y^ Ayjj 

[Id 


Referring to equation (9a), clearly 
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D(x;Ax) = 


exp 


- il (Ax^ + Ay^) 
. AD 




A^D^ 


Tj.(Ax^; x) = exp T- ^ Ax^ + y Ay^)l 
L XD J 

T^(x^; Ax) = (*a ^ 

The other factors must be set approximately equal to unity if the factorization implied by 
equation (9a) is to be possible. The simple interpretation of speckle outlined in this re- 
port requires the factorization implied by equation (9a). Assume 

[Is ^ 

These approximations will be valid provided that the field correlation at the aperture is a 
narrow function of Ax^, Ay^. 

From equation (10) the spatial d^ree of coherence is 

y (x + Ax; x) 


The spatial coherence function y simply involves the Fourier transform of the intensity 
expectation at the object. 

For a tmiform ensemble average of intensity, the correlation of the intensity fluc- 
tiiations stressed by equation (10) is 


= e^ 


[1 i£ (Ax^ + Ay^) 
L AD 


j exp ^ Ax + y Ay)j 


00 

// ya^a'^'a’ V) [I5 <=‘a ?a *y>] <^a “ya 


If 


<®^’^>ya>®a<=V>ya»*a‘’ya 
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X y 


//■ 


rect — rect 


exp 


(x^ Ax + Ay) dx^ dy^ 


Hence, speckle size will be determined by a transform of the shape function of the object. 


= sine 


Ax Lj 
XD 


i sine 


XD 


(19a) 


Referring to equation (13) for a calculation of the area of a typical speckle, 


00 



<-00 


(19b) 


“WALy 

This formula simply expresses the typical area of a speckle, according to one definition 
of area, in terms of the wavelength of light X, the separation D, and the rough-surface 
size by Ly. The actual speckle pattern will appear to consist of irregularly shaped 
speckles of varying size and intensity. 

These results originally derived by Goodman (ref. 8), were derived herein by a sys- 
tematic application of the .analysis developed in this report. The typical extent of a 
speckle corresponds to the first zeros of the sine functions. The power spectral density 
of the intensity fluctuations can be gotten by evaluating the Fourier trai&form of |y|^. 


Speckle Pattern Viewer 

The next sample analysis treats the speckle pattern viewer shown in figure 2. Lens 
L placed two focal lengths from physical aperture stop a creates a unity m^nification 
image of the stop on the television camera lens. The focal length fj^ of lens L is con- 
sidered to be large when compared with the focal length fipy^ television camera 
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lens. The rough surface object is considered to be at a large distance from aperture a. 
Hence, the image of the surface is focused on the vidicon target at approximately f-pyL 
from the camera lens. The impulse response fimction is evaluated in two 

steps. An impulse response fimction is first evaluated from the aperture to the tele- 
vision camera lens. This impulse response function is designated ^L^^V’^TV’^a’^a^' 
Next, an impulse response function is evaluated from the television camera lens to the 
target. This impulse response function is designated hrpy(x,y;x,py,y,py). The impulse 
response function from the aperture a to the target is then 


OO 

hj,^,y;x^,y^) = j I j y j ^yj yTV^^L^^TV’ ^TV’ ^a’ ^a^^TV 


In the analysis that follows, the y -coordinate will be omitted prior to the final results 
since the y-coordinate enters symmetrically with the x-coordinate in the formulas. 

Using reference 7, the impulse response function from aperture to television camera 
lens is 






, 0 0 


/■ 




^ (==TV 


L«L 


dXr 


i <?r (Xt ) 

where Xp^ refers to positions on lens plane L. The pupil function Pl(Xl) = e 
when Xp^ is on the lens, and Pj^(xj^) = 0 when Xj^ is outside the lens, where ^PlCx^) 
is a phase error due to aberrations of the lens L. For this calculation, neglect the ef- 
fect of aberrations so that Pl(xl) = 1 on the lens. Also, assume that the lens is large 
enov^h so that Pl(Xl) = 1 from -<» to «>. Then 


hL(x, 


_ 1 jJTx|y/2XfL 






where 6(x) is the Dirac delta function. 

The impulse response function (ref. 7) from the television lens to the target is 


lVy(x;Xrpy) - — e 

j^fTVL 


jn-x^/Xf, 


TVL „ , -j277X^yX/Xf,j.yL 

iTipy^Xpy;© 
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Again, P^y(x,py) = e on the television lens and Prpy(x,py) = 0 off the tele- 

vision lens. As before, neglect aberrations and set P,j,y(xrpy) = 1 from -<» to °o. 
Because of 6(x,py + x^)/2, hj^ is easily calc\ilated: 


1 J2jrx x/xf^L 


jAf, 


TVL 


Before performing the factorization implied by equation (9a), evaluate 


l^(x + Ax;x^)hjj(x;x^ + Ax^) = e 


A ItVL 


j2?7X^ ^ /^TVL ^^/^TVL 

X 6 © © 


-j2n-x^ Ax/Afrpyj^ j2nx Ax^/Af^yj^ 
X © © 


jn- Ax“/AfT 

The factorization implied by equation (9a) can be performed if e w 1 and if 

i27rx Ax /Af j 

e ^ pa 1 within the aperture or over the interval Ax^^, where the field corre- 

lation is significant. The factorization is 


^ -jTT Ax^/Af^y, -j2;rx Ax/Af„yL 

D(x,Ax) = 1 e ^'^^e 


1 . 2f2 

A ItVL 


Ta(Xa;AS) = e 


“JZttx^ Ax/AfTv^ 


Tc(AXa;x) = e 


j2nx Ax^/Af,j,yL 


Equation (10) then yields the results 
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^ -j27TX Ax/Afrpy^ -jTT Ax /Xfrpyj^ 

y(x + Ax;x)=e ^ e ^ 


r — . — -j27TX Ax/Xfm-^ 

/ <Ea(=^>E^(x^»e 

•^oo 

/ (®^V®a(='a»'*=Si 

rvs 


( 20 ) 


This result is identical with thie expression for y obtained for propagation of light from 
an object of finite size to a detecting plane. The plane separation D is replaced by the 
focal length fr^y^ of the television camera lens. The quantity is the ex- 

pected intensity over the aperture. For example, a circular aperture of radius R with 
approximately uniform intensity in the aperture has an intensity fluctuation correlation 



2ttR 


.Af, 


Ap' 


TVL 


2nR 


Ap 


where Ap 


«/.~2 . ~2 
y Ax + Ay 


The first zero of the correlation that provides an estimate of typical speckle size occxms 
at 

1.22 Afrpyj^ 1. 22 Afipyj^ 

Ap = = 

2R 

a 

where D„ = 2R is the diameter of the aperture. This dimension is identical with the 

A ■ - 

separation of two barely resolved points in the focal plane of a diffraction-limited lens of 
diameter and focal length frj.yL when that lens is used to image a distant object. 


Speckle Interferometer 

The next example treats the laser-speckle interferometer shown in figure 3 (ref. 3). 
This interferometer is but a slight modification of the speckle viewer shown in figure 2 
and described in the previous section. The interferometer is discussed because it ex- 
emplifies an interesting property of speckle: adjacent bright speckles tend to have oppo- 
site phase relative to a reference beam. The single lens L in figure 2 is replaced by 
two lenses as shown in figure 3. A beam splitter is introduced between the two lenses to 
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turn a reference beam toward the target. As shown in figure 3, this arrangement is 
equivalent to placing a point source in the aperture a. Again, the optical parts are as- 
sumed large enough to be considered of infinite lateral extent. As in the speckle viewer, 
the aperture is imaged onto the television camera lens. An analytically equivalent inter- 
ferometer can be formed by placing a point source in the aperture of the speckle viewer 
shown in figure 2 and choosing a lens L with one -half the focal length of the lenses L 
used in the speckle interferometer. The light illuminatii^ the rough surface object and 
the light from the point source are, of course, mutually coherent. The equivalent inter- 
ferometer is shown in figure 4 and is used for the analysis that follows. 

The point source representing the reference beam is represented as E^6(x, y), 
where 6(x,y) is the Dirac delta fvinction.- On the television target in the detecting plane, 
the reference beam is 


= E^h^^;0) = 


Ej, ]7rX Afrpyj^ 

i^TVL 


With respect to both x and y 



jAl, 


exp 


TVL 


in 


Jif, 


TVL 



The formula is the tisual parabolic approximation to a spherical wave. The function 
h|^ was derived in the section on the speckle viewer. K the wave from the rough -surface 
object is designated E^,y), the intensity on the target is 


I(x,y) = 


E I ^ E 

^ . + |E(x,y)|^ - 2 ^ — exp 


x2f2 
A i^VL 


jXf, 


-j_JL_(S2+S^) 


TVL 


Xf, 


E^,y) (21) 


TVL 


The third term in the intensity depends on the relative phase between the reference wave 
and the wave from the rough-sxirface object. Hence, this term is sensitive to a phase 
change in the object wave as a result of object motion. Because adjacent bright speckles 
tend to have opposite phase relative to the reference beam, the speckle size will appear 
to increase when a bright reference beam is present. 

In order to understand these properties of speckle and the speckle interferometer, 
consider the square aperture; 
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Assvuning that the intensity in the square aperture is approximately constant and usii^ 
equation (20), 


+ Ax,y + Ay;x,y) 


-j27TX Ax/XfipY^ -j27ry Ay/Xfrji'yj^ 


-jjT Ax^AfrpYL ^y^/^TVL 

X © © 

From equation (13), the typical area of a speckle can be calculated to be 



Figure 5 shows two typical bright speckles separated by a typical dark speckle. The 
separation along the x-axis between pairs of points, one point in each of the bright 
speckles, varies from to 3Xf,j,yj^/Z. The separation between the same pairs 

along the y-axis varies from 0 to Hence, the average value of the sine func- 

tion product will be negative in equation (22) when comparing points in adjacent bright 
speckles. The phase factors in the first square brackets in equation (22) are a conse- 
quence of Fresnel diffraction. These factors also occur in the reference beam so that 
they do not introduce a phase difference between the reference beam and the object beam. 
However, since the sine function product is negative between adjacent bright speckles, 
adjacent bright speckles differ by a factor e^”^ when compared with the reference beam. 
Adjacent bright speckles therefore have opposite phase when compared with the reference 
beam. 

Referring to equation (21), if one bright speckle results in a positive third term and 
adds intensity to the first two terms, the adjacent bright speckle will tend to contribute a 
negative third term. One bright speckle wiU tend to darken relative to the other. 

Since the speckle interferometer is being used as an example, another property of 
the interference between the speckle pattern and a reference beam will be mentioned. If 
the third term in equation (21) for the intensity contains a time -varying phase factor, 
time averaging may be performed by the detector. Time -average holography, discussed 
for example in reference 7, describes the effect of time averaging on an interference 
pattern. The siun of the first term in equation (21) plus the third term yields an inter- 
ference speckle effect. If the third term has a time-varying phase factor, a time- 
average envelope will be superimposed on the interference speckle. The envelope is a 
function of the amplitude of vibration if the time -varying phase is caused by object 


sme 


AxZ 


Xf. 


sme 


AyZ 


TVL 


Xf, 


( 22 ) 


TVLJ 
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vibration. For sinusoidal vibration, that envelope is the same Bessel function that ap- 
pears in time -average holography. 


CONCLUDING REMARKS 

A general description of the speckle effect requires analytical tools that can handle 
a light source of limited coherence; a reflecting object with arbitrary surface properties 
that manifest themselves in a reflected field of arbitrary spatial correlation, polariza- 
tion, and time variation; an arbitrarily complex intervening optical system that may be 
changing in time; and a detector with arbitrary properties. This report begins with a 
more restricted general analysis that stolid be useful for guiding the design or develop- 
ment of a new speckle technique or instrument. The assumptions are a highly coherent 
light source, the validity of a linear systems approach, possible arbitrary polarization, 
and a time -invariant optical system. This development is presented in equations (1) to 
(8). This analysis is still too general for the simple understanding of, application of, or 
evaluation of existing speckle instrumentation. Since these activities are more common 
in aeronautics applications, the analysis is further developed subject to assumptions to 
give a simple set of equations. The application of equations (5), (9a), (9b), (10), (12a), 
and (13) to (15) more or less in order and the use of geometrical optics should be ad- 
equate to evaluate a speckle instrument or technique. To test this conclusion, a speckle 
viewer , a speckle interferometer , and propagation of light from a rough surface were 
evaluated, . 
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APPENDIX - SYMBOLS 


The variables x, x, and x; subscripted forms of these variables; or increments of 
these variables, when appearing alone as arguments of functions, are understood to 
represent (x, y), and (x,y); subscripted forms of these pairs of variables; or in- 

crements of these pairs of variables. 


s 

D,d 

D^, Ax) 


^TVL 


I(x) 

% 


area of detector 

typical area of speckle 

propagation distances 

fxmction involving possible phase and attenuation factors in detecting plane 

X, y-component of electric field phasor at position (x,y) 

focal length of lens L 

focal lei^h of television camera lens TVL 

probability density function for x-polarized field and y-polarized 
field . Ey at position (x,y) 

probability density function for x-polarized field E^ at position (xj^,yj^) 
and x-polarized field E at position (x2,y2) ^ 

£t 

for deterministically polarized light, probability density function for x- 
polarized field E^ at position (xj,yj) and x-polarized field E^^ at 
position (X2,y2) expressed as function of corresponding y -polarizations 

for deterministically polarized light, probability density function for y- 
polarized field E at position (x< ,y. ) and y-polarized field E at 

yi 11 Y2 

position (X2,y2) expressed as function of corresponding x-polarizations 

probability density function for y-polarized field Ey^ at position (xj^,yj^) 
and y-polarized field E at position (x2,y2) 

for deterministically polarized light, probability density function for y- 
polarized field E at position (x2,y2) and y-polarized field E at 
position (xj^,yj^) expressed as function of corresponding y -polarizations 

impulse response function or point spread function for light polarized in 
X, y-direction 

intensity at (x,y) 

intensity of light in Fraunhofer diffraction; pattern 
intensity of light polarized in x, y-direction 
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j 


imaginary part of complex number 


I 

P(D 

P(x) 

R 

Ta(x^;Ax) 

T^(AXa;x) 

(x,y) 

(2,y) 

^,y) 

<^a>ya^ 

(Xty, Yipy) 

r(x2*,x^) 

rx^y(x2;xi) 

y(x + Ax;x) 
6 

6(x) 

X 

Ap 


lengths of sides of rectangular aperture 
length 

probability density function for intensity I 

pupil fimction (In general, P(x, y) = where (x,y) is a position in 

an apertxjre and P(x,y) = 0 for positions (x,y) outside aperhire; 

^x, y) is a phase function, ) 

radius of circvilar aperture 

function that transforms intensity in limiting aperture to detecting plane 

fxmction that transforms wide-sense stationary function y(Xa + Ax^^jX^) 
from limiting aperture to detecting plane 

position on rough surface 

position on detector 

position in Fraunhofer diffraction pattern of superimposed speckle 
patterns 

position in aperture 
position on thin lens 
position on vidicon 

correlation or spatial coherence between fields at positions (x^,yj^) and 
(x2,y2) when a single polarization is understood to be present 

correlation or spatial coherence between x,y-polarized fields at 
positions (xj,y^) and (x2,y2) 

normalized correlation or spatial coherence function between fields at 
positions (x,y) and (x + Ax,y + Ay); y(x;x) = 1 

relative displacement of overlapping, identical speckle pattern 

Dirac delta function 

wavelength of light 

Ap = Vax^ + Ay^ 
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Subscripts: 
a aperture 

b back part of optical system, or part between limiting aperture and detector 
F Fraunhofer diffraction pattern 

f front part of optical system, or the part between limiting aperture and roTigh sxmface 
t total 
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separated by distance D from an infinite detecting 
plane. 



Beam 

splitter 



Figure 3. - Speckle interferometer. 
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Figure 4. - Equivalent interferometer. 







Figure 5. - Average speckles. 
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